Finite temperature transport in disordered Heisenberg chains 
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Using numerical diagonalization techniques, we explore the effect of local and bond disorder on 
the finite temperature spin and thermal conductivities of the one dimensional anisotropic spin- 1/2 
Heisenberg model. High-temperature results for local disorder show that the dc conductivies are 
finite, apart from the uncorrelated - XY case - where dc transport vanishes. Moreover, at strong 
disorder, we find finite dc conductivities at all temperatures T, except T = 0. The low frequency 
conductivities are characterized by a nonanalytic cusp shape. Similar behavior is found for bond 
disorder. 
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The effect of correlations on localized states in disor- 
dered system^ is a long standing problem 2 , that is attract- 
ing renewed theoretical and experimental interest 3 -^^. 
While it is clear that in a one-dimensional (ID) nonin- 
teracting disordered system all states are localized, the 
introduction of correlations might lead to different pos- 
sible scenaria. First, at zero temperature, T = 0, nu- 
merical results for fermions with repulsive interaction in 
a disordered system reveal that localization persists in 
spite of correlations^, although some types of interac- 
tions might destroy the localized states leading the sys- 
tem to a normal diffusive state or even one with diverging 
low frequency conductivity 10 . Even if the system remains 
localized at T = 0, an arbitrary low temperature could 
delocalize it or a finite critical temperature 4 ^ might be 
needed to drive it to a normal state at high temperatures. 
There are also indications that in the presence of large 
disorder even at high T many body states can appear 
effectively localized^iii. 

To explore this issue, the ID (in general) anisotropic 
spin-1/2 Heisenberg chain is a minimal model that al- 
lows to investigate the interplay of disorder and cor- 
relations on transport. In the XY limit, mapped by 
the Jordan- Wigner transformation to a system of non- 
interacting spinless fermions 12 , it is expected to have all 
single-particle states localized under any amount of local 
or bond disorder, consistent with the Anderson localiza- 
tion phenomenon 1 . The spin (in fcrmionic representation 
equivalent to charge) conductivity as well as the thermal 
d.c. conductivities are expected to vanish at all T, while 
a.c. conductivities are finite but non-trivial^. 

Without disorder, the anisotropic XXZ model is a 
strongly correlated spin system, with nearest neighbor in- 
teraction in the fermion picture. It is integrable using the 
Bethe ansatz method for any value of the anisotropy and 
it is known to show ideal spin/charge (in the easy-plane 



case) and thermal (for any anisotropy) conductivities at 
all temperatures^. 

Besides the theoretical interest of this model quasi- 
1D magnetic compounds have recently been synthesized 
which are described exceedingly well by the ID isotropic 
spin-1/2 Heisenberg model and show unusually high ther- 
mal conductivity due to a magnetic transport mode 
contribution^. Relevant to this work, experiments arc 
underway to study the effect of disorder by non-magnetic 
as well as magnetic impurities. 

In this work, we will use state of the art numerical diag- 
onalization techniques - the exact diagonalization (ED), 
the finite-temperature Lanczos method (FTLM)i^ and 
the microcanonical Lanczos method (MCLM)^ 7 - - to see 
what they can offer on this issue of disorder and corre- 
lations. While we will study the spin and thermal con- 
ductivity of the spin-1/2 anisotropic Heisenberg model, 
the spin conductivity maps directly to that of the charge 
conductivity of a spinless fermion model. 

We first consider the ID anisotropic spin-1/2 Heisen- 
berg model in the presence of a random local magnetic 
field. 
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where S a ,a = x,y,z are spin-1/2 operators, J is the 
magnetic exchange coupling, A the anisotropy parame- 
ter and — W/2 < bi < +W/2 random local fields from 
a uniform distribution. We assume periodic boundary 
conditions, % = Kb = 1 and take J as the unit of energy. 

Our analysis will be based on standard linear response 
theory. The spin conductivity is given by 
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j representing the spin current. The corresponding ther- 
mal conductivity is given by, 

K (uj) = ^-S dte" t <[f(t),f)}>, (3) 

where f3 = 1/T. The energy current j e in an inhomoge- 
neous system can be defined via the dipole operator— 
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whereby hi are local energy operators and n the corre- 
ponding coordinates. Taking into account that locations 
of local field energies are on sites and exchange energies 
on bonds, respectively, one arrives at 
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For vanishing random fields, W = 0, the energy current 
j e commutes with the Hamiltonian for all values of the 
anisotropy A and thus the system is an ideal thermal 
conductor— at all T. It has also been shown that the 
uniform, W — 0, system exhibits ballistic spin transport 
in the XY regime, A < 1, at all temperatures although j 
does not commute with H— . 

Let us start the analysis with the high-T limit, T — > 
oo, which exhibits nontrivial <j{lu) and k(lo), being in 
fact quite generic for all T > 0. The difference between 
spin and thermal transport can be realized already from 
the spin, — J w n o{uj)(Lj = 7t/3to* and corresponding- 
energy, = iT0 2 m' i n frequency moments. Moments can 
be evaluated analytically at T — > oo, e.g. mg = (jj)/L, 
m | = ([H,j][H,j])/L etc. One gets m s = J 2 /8 and 

ml = ^[j2 A 2 +4<fc 2 >], 

= ^[(l + 2A 2 )J 2 + 2{6 2 )], 

m\ = ^(3 + 10A 2 )(fe 2 ) + ^((6 4 )-(& 2 ) 2 ), (6) 

where (b 2 ) = W 2 /12,(b 4 ) = W 4 /80. Eqs.© reveal the 
difference between spin and thermal transport since the 
finite dispersion 5 e = yj m\lm\ of k(lu) is induced solely 
by W > whereas 5 s = y m^/m^ remains finite even for 
W = 0. 

Although the lowest frequency moments can serve as 
a reference, they are insufficient to reveal the most chal- 
lenging lj — > behavior. For the latter we have to rely on 
numerical calculations. The most favorable case for sim- 
ulations on a finite size lattice is the strong disorder limit 



where we expect the localization length £ to be shortest. 
In the following we consider W — 2, where an estimate 
of T = localization length £ exists in the literature^, 
which suggests that £ is less than 10 sites in the cases we 
are studying. 

In Fig. 1 and Fig. 2 we present results for and 
k(uj), respectively. The data for L = 14, with a Hilbcrt 
space dimension of 3432 states in the S z — subsector, 
were obtained by exact (full) diagonalization. The 5- 
peaks at the excitation frequencies are binned in windows 
Slu = 0.01, which also gives the frequency resolution of 
the spectra. There is an average over N r = 10 random 
field configurations. For L = 16—24 the MCLMii is used, 
particularly suitable for high T(^> J), with typically 2000 
Lanczos steps and random-configuration sampling N r > 
100. In the same plots, we show the noninteracting case 
for L = 1000 and averaging N r = 1000, where we expect 
the dc conductivities a ( ic — cr(u> — > 0) and Kdc — k(cj — > 
0) to vanish. All the spectra are normalized to a unit 
integral. 




FIG. 1: Dynamical spin conductivity <j(uj) at T — > oo for local 
disorder W — 2 and various A (curves normalized to unity), 
evaluated via ED (L = 14) and for A = 1 also via MCLM 
(L = 24). 
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FIG. 2: Dynamical thermal conductivity k(lu) at T — > oo for 
W = 2 and various A (curves normalized to unity) , evaluated 
via ED (L = 14) and MCLM (L = 24) 
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Data for various sizes L indicate a convergence of the 
finite-size effects at moderate sizes L > 16 (at least for 
A > 0.5) for the chosen rather strong disorder W — 2. 
In particular, we show in Fig. 1 comparison of L = 14 
and L = 24 results for A = 1. In Fig. 2, for the same 
parameters, the curves are also nearly indistinguishable 
for L between 14 and 24 sites. 

These results clearly reveal that apart from the XY 
limit, A = 0, the system is conducting, i.e. both spin 
Ode and thermal Kdc d.c. conductivities are finite. Nev- 
ertheless, due to the large disorder W the dynamics is 
non Drude-like, since the maximum of u{lo) as well as 
of k(uj) appears at a finite lu* > 0, in analogy with 
the localization at A = 0. Hence, at A > and large 
W we are dealing with pseudo-localized dynamic a 6 ' 11 . 
Another novel feature of this regime appears to be a 
generic (nonanalytic) cusp-like behavior at low frequen- 
cies, cr(u)) ~ &dc + ot\u}\, k(lo) ~ Kdc + tM, for which so 
far we cannot offer an analysis. It might be attributed 
to long-time tail effects although, in such a case, the low 
frequency drop of the conductivity was found to be only 
a few percent^ and not by an order of magnitude as 
in our case. Such a frequency dependence is strongly 
reminiscent of the behavior in strongly disordered 2D 
system as has been analyzed theoretically and observed 
experimentally 21 . 

Apart from a qualitative similarity between cr(u>) and 
k(u>) in Figs. 1,2 there are also some differences. er(w) 
is more sensitive to A, as it is already evident from the 
moments, Eqs.(|6]), and corresponding 5 s . They originate 
from the fact that even at W = cr(u> — > 0) changes 
qualitatively at A = 1, not being the case for k(lo — > 0) 14 . 

As our numerical simulations indicate, a similar qual- 
itative behavior persists by decreasing the disorder to 
W = 1.0 (not shown). With decreasing W the pseudo- 
localized form gives way to a more Drude-like form with 
lu* — > and strongly increased ddc, Kdc- However, reduc- 
ing the disorder further we are running to long localiza- 
tion lengths (for A = 0) and in general less controllable 
finite size effects preventing reliable conclusions. 

The next issue is the temperature dependence of the 
dynamical (in particular dc) conductivity and the even- 
tual existence of a critical temperature T c below which 
the system becomes insulating 4 ^. To study this question 
we employed the ED method for L = 14 (using N r = 10) 
and the FTLM for L = 16 - 20 with 200 - 400 Lanc- 
zos steps for high frequency resolution and N r ~ 100. 
The results are qualitatively similar whereby the FTLM, 
properly interpolating between the T — (ground state) 
Lanczos method and T > behavior, is more reliable for 
small T < 0.5 due to larger L and more dense low energy 
spectra. Results for a{u>) and k(uj) in the isotropic case 
A = 1 and at fixed W = 2 are shown in Figs. 3, 4 for 
various T = — 2 and L = 20 (being essentially equal to 
the results obtained for L = 16). The data again indicate 
that Ode and Kdc remain finite at all T > vanishing only 
at T = 0. A rather abrupt drop of Ode appears at T ~ 0.1 
which is however in the range of finite-size temperature 



T fs (for available L = 20) below which the FTLM results 
are not to be trusted 1 ^. These data suggest a zero criti- 
cal temperature of localization-delocalization transition, 
although of course we cannot exclude an exponentially 
small one, which is beyond the reach of actual numerical 
simulations. 




FIG. 3: Spin conductivity a(co) for A = 1 and W = 2 for 
various T. 
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FIG. 4: Thermal conductivity k(lo) for A = 1 and W = 2 for 
various T. 

In connection with existing ID magnetic compounds 15 
more relevant appears to be the spin- 1/2 (anisotropic) 
Hamiltonian with bond disorder, i.e. disorder in ex- 
change couplings, 

H = J2 Ji,i+i(SfSf +1 + SfSf +1 + ASf Sf +1 ), (7) 
i 

where J/,;+i = J(l — and we assume —W/2 < 

Sij+i < W/2 uniformly distributed random numbers. 
Such disorder can be induced, e.g., by coupling to static 
lattice displacements^. The local spin current is now 
jy+i = Ji,i+i(Sf Sf +1 - SfSf +1 ) while the energy cur- 
rent is given by 

j e = X>-i,^y+iif • (8) 
i 

An open question is whether a ID spin chain with bond 
disorder 10 behaves qualitatively different to the site dis- 
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order discussed above. Our results indicate that it is not 
the case. 
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FIG. 5: T — > oo results for a{u>) and for A = 1 and different 
disorder W (curves are normalized). 
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FIG. 6: T -> oo results n(u) for A = 1 and different W 
(curves are normalized). 

In Figs. 5, 6 we present T — ► oo results for u(ui) as well 
as k(uj) for A = 1 and different bond disorder strengths 
W = 0.5, 1, 1.5. Results were obtained using the MCLM 
method on L = 20 sites. Results for larger W — 1, 1.5 are 
well converged with size and clearly indicate that we are 



again dealing with finite dc limits Ode > and Kd c > 0. 
With respect to the site disorder case in Figs. 1, 2 there 
are similarities but also differences: a) for bond disorder 
we are restricted to W < 2 to have a meaningful model 
without a possibility of a broken bond, b) the pseudo- 
localization is less pronounced at least for a(uj) and shows 
up only closer to W = 2, e.g. for k(uj) at W — 1.5, c) k(oj) 
in Fig. 6 reveals a quite abrupt crossover with disorder 
strength, from a Drude-like response (at W — 0.5) to a 
localized-like one with ui* > at W — 1.0, d) at least for 
cr(u>) two energy scales are evident in Fig. 5 which are 
not present in the random-field case. 

In conclusion, our results of numerical simulations on 
the interplay of disorder and correlations in the spin and 
thermal transport within Heisenberg spin chains can be 
summarized by the following scenario: (a) finite random- 
field disorder W > induces localization and vanish- 
ing dc transport at any T in the XY limit, correspond- 
ing to noninteracting fermionsi, and as well generally at 
T = 0&2 (for A > considered here); (b) apart from 
the latter two limits the system appears to behave as a 
normal conductor with finite a^c > 0, Kdc > both for 
various A > and T > 0; in particular, we do not find 
any evidence for a phase transition by varying T or W; 
(c) dynamical transport (at least for larger disorder) re- 
veals a generic cusp- like nonanalytic behavior for ui, anal- 
ogous to long-time tails in classical dynamical systems in 
low-dimensional- 3 or 2D strongly disordered systems; (d) 
with increasing disorder the system reveals a crossover 
from the Drude-like to a pseudo-localized dynamics with 
very low dc <Jd c , Kdcr*j an d (e) similar conclusions seem 
to hold for the bond disorder. 

Clearly, several caveats are in order. The considered 
cases mostly correspond to substantial disorder, where 
the finite size effects are well under control and results 
converged within available L, at least for T > Tf s and not 
too small A > 0. Also, numerical results cannot exclude 
the localization on a very long scale £ 3> L although we 
do not find any signature of such a development. 

This work was supported by the FP6-032980-2 NOV- 
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